Decoherence of collective atomic spin states due to inhomogeneous coupling 
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We investigate the decoherence of a superposition of symmetric collective internal states of an 
atomic ensemble due to inhomogeneous coupling to external control fields. For asymptotically large 
system, we find the characteristic decoherence rate scales as y/N with N being the total number 
of atoms. Our results shed new light on attempts for quantum information processing and storage 
with atomic ensembles. 
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Coherent quantum information encoding and process- 
ing has recently emerged as a major goal for the physics 
community. Despite the seemingly insurmountable diffi- 
culties, a rich variety of implementations are being pur- 
sued in laboratories across the globe. Among the early 
success is the apparent ability in simulating quantum op- 
erations with liquid based NMR jl], . Recent theoretical 
efforts indicate, however, that in the pseudo-pure state 
approach using NMR, quantum entanglement, a key ele- 
ment for powerful quantum information processing, was 
in fact no present || [| . Room temperature NMR tech- 
nique, is therefore limited and can not be explored fully 
to benefit from an exponentially large Hilbert space of 
only a polynominally scaled resources and controls ||. 
Nevertheless, early NMR based experiments have pro- 
vided useful insight into the operations of genuine quan- 
tum computers H. 

Over the last few years, using symmetric collective in- 
ternal states of an atomic ensemble, has attracted much 
attention [Q, 0, ^|, [h], |ll). The pros and cons of such 
an approach assisted with cavity photon-atom interac- 
tion was recently discussed by Fleischhauer et al. [l0| . 
In normal cavity QED based quantum computing imple- 
mentations, atomic qubits are entangled and logic opera- 
tions performed through their interaction with the com- 
mon cavity photon quantum field. To maintain quantum 
coherence, it is important to reach the so-called strong 
coupling regime, when the single-photon coherent cou- 
pling go 3> 7j k j the atomic and cavity dissipation (de- 
coherence) rates respectively ]l2| . The symmetric collec- 
tive internal states can reach the strong coupling regime 
without requiring a high finesse cavity as go oc V^/V, with 
N the total number of atoms jy^, |lj, |l5) . When imple- 
mented with protocols insensitive to individual atomic 
dissipation/decoherence rate as in the Dark state based 
adiabatic transfer protocols [fflpj, |l6||, one can apparently 
gain as upper hand over systems based on single atoms 
inside a cavity ||, |o). This is in fact, not quite surpris- 
ing, earlier cavity QED experiments have relied on the 
enhanced dipole interaction of a collection of many atoms 
p3| , [l4|, |l5). In free space, the phenomena of superfluores- 
ccncc or super-radiance [TtJ constitutes another ex- 
ample of collective state dynamics. Recent experimental 
success clearly demonstrates the power of such an atomic 
ensemble based system for entangling macroscopic ob- 
jects llS, n9[. Several new ideas have raised further ex- 



pectation of exciting developments to come [£C| |2lJ . Nev- 
ertheless, all ensemble based systems suffer from the re- 
duced size of computational Hilbert space. In this case 
of symmetric collective internal states, the space used for 
quantum information, is much less than the V T 



for N two level atoms [22 . In view of the recent experi- 
mental success in storage and recovery of light coherence 
in atomic gases |23, 24 , a related question to address is 



the sensitivity to errors when collective spin states are as 
quantum memories. 

In this paper, we investigate the decoherence for a su- 
perposition of symmetric internal states of an atomic gas 
due to its inhomogeneous coupling with external control 
fields. Our study is motivated by the simple observation 
that the symmetric states of an atomic ensemble spans 
the computation space only if atoms can be manipulated 
cooperatively, namely, the coupling of both the exter- 
nal manipulating field and the environment surrounding 
the atomic ensemble should be homogeneous such that 
the collective motion of the atomic ensemble can be de- 
scribed by the collective quasi-spin operators. In essence 
the effect of different spatial positions for atoms 1, 2, • • •, 
and N, is ignored or absorbed into each single spin oper- 
ators. In reality, an optically thick atomic ensemble suf- 
fers from inhomogeneous coupling to both classical and 
quantum light fields, i.e. the coupling strength is po- 
sition dependent. Such a situation arises naturally for 
trapped ions due to its center of mass motion. In this 
case, it is well known that the loss of quantum coher- 
ence for a superposition of internal state occurs. In this 
study, we refer such decoherence effect as inhomogeneous 
decoherence. We focus on introducing our technique and 
study the simplest example of superpositions of collective 
atomic Dicke states in this paper. The consideration of 
Dark-state, polariton approach based proposals (l(J will 
be given in the future, as significant complications arise 
when the quantum cavity field is included. 

Our model constitutes of an ensemble of two-level 
atoms described by Hamiltonian 



(*) , I 



(k) (fc) 



h.c.% 



(1) 



k=l 



where the cr's are the standard Pauli matrices (Ti = 1), 
and the different local coupling g^ (for the fc-th atom) 
may be due to a cavity mode profile as common in tightly 



2 



focused cavities or when atomic motional wave packet becomes the focus of our study with 



is insufficiently localized |25|. u>i^ = e y a' -ul, is the 
difference between atomic energy £„ and the external 
near resonant laser frequency ojl- For convenience, we 
further abstract Eq. ([j]) into the compact form H = 

EfcLi B (k) ■ v (k \ with real parameters B^K 

The symmetric collective spin space V s of dimension 
2,7 + 1 < 2 N {J = N/2) is spanned by the collective 
angular momentum states {| J, M), M = — J, ■ ■ ■ , J— 1, 7} 
of J i* = J2iLi a ^ / 2 satisfying [J M , J v ] = ie^J^ and 
J% + Jy + J% = J 2 = J (J + 1). e M „c IS tnc symmetric 
permutation tensor. The \ J,M) space can be generated 
from the ladder operator J± = J x ±iJ y according to [p6[, 



.(A) 



\J,M) 



(J-M)! 



rJ+Af 



(J + M)!(2J) 



(2) 



except we note that an arbitrary unimodular pha- 
sor can be self-consistently included with J± = 



E£Li e ±l ^a^/2 and | J, - J) = | • • • , |). 

For any realistic system, an inhomogeneous distribu- 
tion of the parameter BV> makes it impossible to con- 
strain the system dynamics within the subspace V s . To 
facility further discussion denote H = Ho + Hi with 
H a = EfcLi B ■ ctW and Hi = Ef=i b (k) ■ ff {k \ where 
£?( fc ) = B+b^ with B = J2 k B [k) /N. H constitutes the 
intended coupling between the symmetric collective spin 
states, while Hi represents a source of inhomogeneous de- 
coherence. It causes decoherence as it provides a direct 
coupling from the subspace V s to its complement V° in 
V T . A quantitative measure for the unwanted coupling 
Hi is in terms of the leakage parameter. Suppose initially 
the system is prepared in a superposition of collective 
spin states |^(0)) S V s . The intended dynamics gov- 
ltHo = nf=i e- ltS -° ik) leads to the 



erned by Uo(t) = e 



resultant state \4>{t))o = Uo{t)\<p(Q)) , still within the same 
subspace. The actual final state is \4>(t)) — U(t)\<fr(0)) 



with U(t) = fl 



N 

k=l • 



generally will span more 



than V s . The leakage can therefore be defined as 

e =i-\(Mt)m))\ i 



(3) 



£ = corresponds to no leakage, while £ — > 1 indicates a 
complete loss of the system coherence and population. 

Denote \4>(0)) € V s as a normalized state expanded in 
terms of | J, M) , 



10(0)) 



c M \J,M), 



(4) 



M<iV,orM~iV 

the overlap 

KMt)\m\ 2 =\{<f>(o)\ukw(t)\m)\ 2 

= J2J2 C M' C MO M 'M(t)<l, (5) 
M M' 



(J,M'\U^t)U(t)\J,M) 

N 

(J,M'\l[0to\J,M), 

k=l 



(6) 



0(A) = R( k ) +if( k ) .^ k \ and 

i? (fe) = cos Bt cos B (k h + (n ■ n (fe) ) sin£tsin5 (fe) t, 
7 (fc) = h sin Bt cos B (k) t + n (k) cos Bt sin B (k k 

+ (n x h (k) ) sin B sin B (k) t. (7) 

We have defined n = B/S and riW = 

The evaluation of Eq. (^) is difficult as state | J, M) 
involves a symmetric permutation of all atoms so that 
the FJ fc factor can not be pulled outside the inner prod- 
uct. Furthermore, IlfcLi expands into 2 N separate 
terms, involving asymmetric products of (?( fe ) of upto 
powers of N. A similar product structure was found to 
be responsible for decoherence in quantum measurement 
models pfl] , where the decoherence factor (the overlaps of 
the final states of detector or a environment) suppresses 
the off-diagonal element of its reduced density matrix. 
In mathematical terms, for a factorized state |/) = 

nf=i l/ (fe) >, the overlap integral (/| Uk=i W$ M \f) be- 
comes n^i(/^I^M'Afl/^)> wn i cn approaches zero 



in the limit of 



N as each factor 



(f(ty\wffi M \fW) has a norm less than unity. To make 



macroscopic 

. m'm\J 

a similar argument for the present problem, we need to 
find an expression such that the collective state | J, M) 
becomes factorized. Since we are interested in obtaining 
the asymptotically valid results in the limit of large N, 
a short time approximation (small t) can not be simply 
adopted. Following early discussions on atomic coherent 
states [B8L E91 , we introduce 



a phase coherent state, that can be expanded according 
to the number of excitations 



1 



2 N/2 



l + e ie J+ + ••■ + — Jl 



J- e KJ+M)6 



M=-J 



JM 



| J, -J) 

(9) 



where Njm = V( J + M )-(J - M)\2 N /{2J)\. The in- 
verse transformation gives 



|J,M) = 



AT. 



JM 



e -i(J+M)0 (1Q) 



2vr Jo 

which helps to evaluate Eq. @ as Om'm = 
NjmNjm'Om 1 m with the reduced overlap 
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FIG. 1: M and J dependence of / for J = 500 and B z 
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e -i(J+M)e e i(J+M')B' 



J\G {k) {6,6') (11) 



k=l 



in a simple factorized form and 



G (k) 



1 



-CO 



,UI(i 



■)OW(l 



-)U>* 



|G«| < las both (l+e i9 ^)l 4)fc/\/2 and (l+e ie 4^)1 i 
)fc/\/2 are normalized. This points to a strong physi- 
cal argument against rapid decoherence of collective spin 
state qubits. The question to answer is now clearly how 
does Omm approach due to inhomogeneous coupling. 
If the coupling coefficients and u>a were constants 
(independent of atom label k), Omm 1 = Smm>- 

We investigate the above question for several model 
cases. First, we look at inhomogeneous broadening when 



(fc) _ B (k) _ 



and B 



satisfies a normal distribution 



(with respect to k) with mean B = B z z — (B z )z, and 
variance a 2 . We find Omm' (t) oc &mm' with the coeffi- 
cient being a constant unity for \M\ = J but decays with 
a time constant T\/ 2 oc l/(\^Na z ) for \M\ < J . Define 
T X j2 = l/(/cr 2 ), we find / is essentially independent of 
a z for a z £ [10~ 7 , 10 _1 ]-B Z . It contains an apparent de- 
pendence on J 2 — M 2 as shown in Fig. |l| for a given J 
and B z . The J dependence (for M = 0) is also shown in 
the same figure. Based on our extensive numerical study, 
we find to a high level of accuracy 



T x/2 {J,M,a z 



1 



no z VJyJl- M 2 /J 2 ' 



(12) 



with k (~ 1.2), essentially independent of B z for B z E 

[icr 2 ,l0 2 ]. 

This result is to be expected based on the collapse and 
revival of a quantum wave packet fl30|] , since each indi- 
vidual atom collapses with a time constant oc 1/<7 Z , the 
collective states of an Guassian ensemble should collapse 
with a time constant oc l/(^f~N a z ) as the net variance 
simply adds. This is indeed what we find for M — or in 
general for \M\ <C J. Equation ( |l2"| ) also indicates that 
significantly reduced decoherence does occur in this case 
for \M\ ~ J, a regime where collective spin states are 
mostly useful § g, ^ |(J In fact, for a single 



FIG. 2: Periodic behavior for \OMKi(t)\ ■ Solid line denotes 
M — 4 and N = 10 with respect to the lower time axis, while 
dashed line denotes M = 2 and N = 2; — k is taken 

for iV = 10 to assure the appearance of revival. For N = 2 
revival occurs for arbitrary random values of B z . 
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FIG. 3: The M dependence of / for J = 200 and B r = 10. 
The smooth curve is a fit given by -4.06483 x 10" 7 |M| 3 + 
2.03393" 4 M 2 + 0.00697|A/| + 10.62188. 



qubit quantum memory involving the two state superpo- 
sition of M = — J and — J + 1, the decoherence rate is 
just that of a single atom ||. For small values of N, when 

(k) 

ratios of different coupling strength B z match ratios of 
integers, we indeed were able to find the expected revival 
as shown in Fig. |^. This of course will not happen for 
an ensemble with a macroscopic N. 

Next we consider the case of inhomogeneous Rabi cou- 

(k) 

pling with B x y being Gaussian distributions with mean 



B 



B„ = B r and variance a 2 = a 2 



t 2 and B { z k) = 0. 



Similar to the previous case, we find the diagonal term 
OmmW (including M — ±J) decays with a time con- 
stant T1/2 = 1/ f<7 r . The J dependence of T1/2 is in fact 
almost identical, i.e. Jm=o — «iJ 1/2 , with k\ w 0.76 
when B r = 10. The M dependence, on the other hand 
is more complicated as shown in Fig. |. Obviously / 
does not depend on M linearly as now Omm(^)| seems 
to decay faster for larger values of \M\. 

The off-diagonal element Om/m'W grows to signifi- 
cant nonzero values, as shown by the typical sampling 
of \OMM'{t)\ 2 in Fig. ^ when N is not too large. Over- 
all, we find the dependence on the random number sam- 
pling is strong only when M — M' — ±1, so we focus 



on M - M' = ±2 here. Define T„ 



as the time for 



\Omm' (t) I to reach its first maximum and O max the value 
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FIG. 4: Typical sampling of O mm > (t) for J = 100, M = J— I, 
and M' = M-l (solid line) and M-2 (dashed line). Different 
figures correspond to different random number sets. 

of the maximum. We find that similar to the diagonals, 
T max = 1/ fcr r , with / a function of J, M, M', and B r , al- 
though O max seems to be largely independent of ay. To 
study the J dependence of / and O max , we consider the 
limiting case of \M\ ~ J when collective states are usual 
proposed to work. The result / oc J 1 / 2 is once again 
as expected. In this case, we also find quite accurately 

^max J 



To summarize, we find within our model, the appar- 
ent decoherence or dissipation rate for superpositions of 
collective spin states scales as \/N. This evidence clearly 
demonstrates that asymptotically there is no advantage 
of using collective spin states for quantum information 
processing. The V7V enhanced coherent dynamics is sim- 
ply being compensated by the y/~N enhanced decoherence 
when inhomogeneous coupling arises. 

Finally, we note that our result also applies to the 
case of entangled states between the collective spins of 
two separate ensembles. For instance, for two ensembles 
A and B, a state Y.m a m b c m a ,m b \Ja, M a ) a \J b ,M b )b 
can always be expressed as coherent superposition of 
the total angular momentum basis J = J a + Jb, i- e - 
into collective basis \Ja, Jb) J,M — Ma + Mb) with 
J = J a + Jb = (N A + N B )/2. 
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